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I. Introduction, work motivation and summary of contributions 

In recent years, considerable attention has been given to the issue of joint transmitter and receiver optimization 
in wireless communication systems. Indeed, many studies have shown that coupling an optimized receiver structure 
with the use of judiciously designed signals at the transmitter can result in substantial performance improvements 
with respect to the case in which transmit signal optimization is not taken into account. 

With reference to direct-sequence/code-division-multiple-access (DS/CDMA) systems, the leading air-interface 
technology of current third-generation wireless data networks, transmitter optimization essentially amounts to the 
problem of spreading code allocation among users. This has been for years and continues to be a very active 
research area. Among the pioneering contributions in this field we find [1], where the sum-capacity maximizing 
sequences for an equal power synchronous DS/CDMA system are derived. The results of this study are then 
reconsidered and generaUzed in several other papers. In reference [2], as an example, the equal power constraint 
of [1] is removed, and optimal sequences for sum-capacity maximization in a synchronous system with arbitrary 
power profile are derived. The paper [3] considers instead the problem of sum-power minimization under pre- 
assigned quaUty of service (QoS) constraints, and shows that, when the sum-power minimizing signatures are 
used, the multiuser Unear minimum mean square error (MMSE) receivers reduce to conventional matched filters. 
In [4]-[7], instead, iterative algorithms that converge to these optimal spreading sequences are given, and it is also 
shown that maximization of the sum capacity is equivalent to minimization of the total squared correlation and 
to the minimization of the global (sum) mean square error (MSE). Similar results are also reported in [8]. The 
papers [9]-[ll], instead, generalize part of these results to the case in which the multiuser system is affected by 
multipath fading. While the above cited papers consider the problem of spreading code design aimed at optimizing 
a global performance measure, such as the sum-capacity, the sum-power, and the global MSE, the focus of this 
paper is on competitive (i.e., non-cooperative) joint spreading code and receiver optimization, wherein each user 
selfishly chooses its own spreading code aimed at optimization of an individual performance measure, such as, 
for instance, the MSE incurred by that user. Game theory [12], a mathematical theory suited to describe the 
interactions among entities with contrasting interests, has emerged in the last decade as a natural tool to study 
such distributed resource allocation algorithms [13], [14], since it provides a natural framework for the design and 
analysis of non-cooperative behavior. Indeed, the selfish users' interaction in a wireless network can be modeled 
as a non-cooperative game in which the users' terminals are the players of the game tuning parameters such as 
transmit power and adopted spreading code in order to compete for network resources and to maximize a given 
utiUty measure. Any action taken by a user, i.e. the choice of the waveform to transmit and of the radiated power, 
modifies the multiuser interference scenario in the network, and thus affects the performance of other users in 
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the network. Game theory is the natural tool for studying this kind of interaction. 

A game-theoretic framework for non-cooperative energy efficiency maximization has been widely applied in 
the recent past to design resource allocation policies for DS/CDMA systems [15]-[17], ultrawideband (UWB) 
systems [18], and multiuser multiple-input multiple-output (MIMO) systems [19]. In particular, the paper [15] 
focuses on the problem of power control for non-cooperative energy efficiency maximization in a CDMA system. 
It is shown that the considered problem admits a unique Nash equilibrium (NE), and, also, pricing is used to 
force the NE point to be Pareto-efficient. The results of [15] are extended in [18] to the case of a multiuser UWB 
system subject to a multipath fading channel, and in [16] to the case in which energy efficiency is maximized in 
a multiuser CDMA system with respect to both transmit power and choice of the uplink linear receiver. Putting 
together the results of [16] with those regarding the issue of spreading code optimization [4]-[7], in [17] the 
problem of energy efficiency maximization with respect to the choice of the transmit power, uplink linear receiver 
and transmit spreading code is addressed; it is shown here that the considered problem admits a unique NE point, 
that, remarkably, is also Pareto-efficient under the assumption that the number of active users does not exceed 
the processing gain. Similar results are also reported in [19] with reference to a multiuser MIMO system, while 
the problem of adaptive implementation of the games developed in [17] is considered in the study [20]. 

All of these works, however, consider, for the sake of simplicity, a real signal model. Otherwise stated, I and 
Q components are not taken into account; while at a first look it might seem that the extension of the reported 
results to the complex case is straightforward, this is not always the case, especially when considering issues 
regarding the problem of spreading code optimization, such as for instance those reported in [4], [6]. 

On the other hand, it is known that if modulations such as binary phase shift keying (BPSK), multi- 
level amplitude-shift-keying (ASK), Gaussian minimum-shift-keying (GMSK), and offset quadrature-amplitude- 
modulation (QAM) are used, receiver performance can be improved through the use of widely-linear (WL) 
reception structures, i.e. receivers that separately process the data and their conjugates [21]-[25]. Indeed, when 
the modulations cited above are employed, the baseband equivalent of the data is an improper complex random 
processQ. 

This paper is thus devoted to the problem of non-cooperative resource allocation in a multiuser DS/CDMA 
wireless data network coupled with WL filtering. In particular, we consider the non-cooperative maximum-SINR 
game for spreading code allocation and WL receiver choice. A distributed iterative algorithm is proposed for 
signature and receiver update. This algorithm is shown to be convergent and to admit one unique stable fixed 

'a complex random process z{t) is said to be proper if its pseudoautocorrelation function Rz{t, u) = E {z{t)z{u)'\ {E{-} denotes the 
statistical expectation) is zero Vt,it while is said to be improper when Rz{t,u) is non-zero. 
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point. Several properties of this algorithm are also shown, mainly generalizing to the case of WL filtering and 
complex signal model many results scattered in references [2], [4]-[7]. A non-cooperative maximum energy- 
efficiency game is then considered, showing that also in this case the use of WL filtering brings a substantial 
performance improvement with respect to the case in which linear detection is used at the receiver side. The 
analysis of this game is also carried out through large system analysis (LSA) arguments. 

The rest of this paper is organized as follows. The next section contains the system model and a brief review 
of linear and WL detection rules. Section III is devoted to the problem of optimal spreading code allocation 
for non-cooperative SINR maximization: an iterative algorithm is proposed which is shown to admit one unique 
stable fixed point, corresponding to the global minimizer of the total mean square error. Section IV deals with the 
problem of non-cooperative energy efficiency maximization, assuming that also the transmit power is among the 
parameters to be tuned. Section V contains the numerical results, while concluding remarks are given in Section 
VI. 



II. System model and multiuser detection structure 

Consider the uplink of a K-user synchronous, single-cell, direct-sequence (DS)/CDMA network with processing 
gain N and subject to fiat fading. In each symbol interval, the received signal is chip-matched filtered and sampled 
at the chip-rate, yielding the A^-dimensional received vector 

K 

r = YjVPkhke^'^''bkSk + n , (1) 

k=l 

where pk is the transmit power of the A;-th user, bk G {—1, 1} is the information symbol of the A;-th user, h^e^^'' 
is the complex channel gain, and is the spreading code of the k-th user. In the sequel we assume that Sk 
is constrained to have unit norm. Finally, n is the thermal noise, assumed to be a zero-mean white Gaussian 
random process with covariance matrix IMqIn, with the identity matrix of order A^. 
Now, a linear detection structure decides on the symbol according to the rule 

h = sign{^[dj^r]} , (2) 

with (•)^ denoting conjugate -transpose, bk the estimate of bk, the A^-dimensional vector representing the 
receive filter for the user k and 3f?[-] denoting real part. The signal-to-interference-plus-noise ratio (SINR) 7^ 
corresponding to the decision rule Q can be written as 

^ Pkhl\d^su? 
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It is well-known that the SINR-maximizing receiver is the linear MMSE receiver, namely 



(4) 



K 



with M = E {rr^} = '"^^pkhls^sjf + IMqIn the covariance matrix of the data. 

k=l 

On the other hand, a WL detection structure can be also devised, i.e. the following detection rule can be 
considered: 

h = sign {5? [df^r + d^^^r*] } , (5) 

with (•)* denoting complex conjugate and where di ^ and d2,fc are two A^-dimensional vectors. Defining the 
augmented vectors 



K 



y^ihibiSi^a + na , and d^.a 



1=1 



d2,k 



(6) 



with 



n 



and Sk^a = 



Ske- 



(7) 



the SINR corresponding to the decision rule ([Sjl is expressed as 



7fc 



2^^o\\dk,af + Y,'^Pihi\dk,a^^,a\^ ' 

iytk 

The SINR maximizing receiver is now written as 



(8) 



dk,a = V^PkhkMa ^Skn, 



(9) 



K 



with Ma = E {vaTa^ = ^ 2pkhlsk,aSk,a + '^o^2N the covariance matrix of the vector r^. It is easily shown 

k=l 

that 



Ma 



M M' 

M'* M* 



(10) 



K 



with M' = '^^Pkhle^'^^'' SksJ^ the pseudo-covariance matrix of the data vector r. The MSE achieved by the 



k=l 



k-th user, is now defined as 



MSEk = E{\bk-dj^^ara\^} , 



(11) 
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If receiver ^ is used, it is easily shown that the k-th user's MSE is written as 

MSEk^E{\bk-d^,ra\^} = , (12) 

wherein 7^ is the k-th user's SINR when receiver Q is used. 

III. MiNIMUM-MSE NON-COOPERATIVE JOINT SPREADING CODE AND WL RECEIVER OPTIMIZATION 

Equipped with the above data model, in this section we consider the issue of non-cooperative joint spreading 
code and WL receiver optimization; the objective function here considered is the individual SINR and/or the 
individual MSE, since, as we will see, the problems of SINR-maximization and of MSE-minimization will turn 
out to be equivalent. 

We thus start by considering the game 

min MSEfc , A; = 1, . . . , - 1 , (13) 

wherein s^, di ^ and d2,fc are A^-dimensional complex vectors, with constrained to have unit norm. 

To analyze this game, we first note that a joint minimization of the MSE with respect to the receive vectors 
and to the transmit spreading code cannot be done in closed form, so we resort to the usual iterative approach 
wherein MSE is minimized sequentially with respect to the receive vectors and to the spreading codes, and, after 
convergence has occurred, make an a posteriori check that the global minimum of the cost function has been 
achieved. We have already noted that the WL receiver for the k-th user minimizing the k-th user's MSE in ([TTI l 
is given by the MMSE receiver in 

Minimization of ([T2l ) with respect to s^ ^, subject to the constraint that ||sfc,a|P = 1> is straightforward, and 
gives the solution 

Sk,a = II , II • (14) 

ll"fc,a|| 

Summing up, an iterative non-cooperative algorithm for MSE-minimization is the following: 



dk,a = V^hkMa ^Sk,a ^k = 1, . . . ,K , 



sfc,a = p^ yk = i,...,K. 



(15) 



Iteration (fTSl ) resembles the well-known MMSE signature update algorithm 

dk = ^hke^'P'M-^Sk yk = l,...,K, 

^k = ST yk = i,...,K, 



(16) 
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which was first proposed in [4] and has been thoroughly analyzed in a number of subsequent studies such as [5]- 
[7]. As already noted, these results apply to the case of a real signal model with linear detection rule; iterations 
([TS] ). instead, come from the non-cooperative MSE minimization assuming a complex signal model and WL 
processing at the receiver. The sequel of this paper is thus devoted to the analysis of the properties of the steady- 
state solution of (fTSl) and to a performance analysis in comparison with the fixed point of iterations ([T6l) . Our 
algorithm starts with K unit norm signatures [si ,j(0) . . . sj^ ^(O)] that we collect into the matrix <S'a(0). At iteration 
n the algorithm replaces each signature with its normalized MMSE filter, thus building the matrix Sain). Notice 
that each iteration is composed of K intermediate steps. At the A;-th intermediate step in iteration n the first k — 1 
signatures have been updated to yield the matrix = [si^a(n) . . . S(fc-i),a('^)) ^A:,a("- — 1) • • • SK^ain — l)]. 

Therefore the interference plus noise correlation matrix employed to update the signature s^^a is 

^k{n) = J2j<k '2pkhk^Sj^a{n)sl^{n) + Y.j>k 2pkhk^Sj,a{n - - 1) + 2N0I2N ■ ^^^^ 

In the following we prove important theoretical results for the algorithm ([TSl l. 

A. Proof of convergence of iterations f li5l ) 

We start by showing that iterations ([T5] ) always converge, i.e. that an NE exists for the considered game ( fT3] ). 
Recalling that a WL receiver operates according to the decision rule (IS), and that its performance is ruled by 
the real part of the cross-correlation matrix of the adopted complex spreading codes, we introduce the WL total 
weighted square Correlation (WL-TWSC) as follows: 

K K K K 

WL-TWSC = ^Kp,/i2/j2(3fJ(sf = Y,Y.PiP^^i^M-^^'-\'' ■ (1^) 

i=l j = l 1=1 j = l 

Since the WL-TWSC is lower bounded for any value of K, the convergence of iterations ( fTSl ) can be easily 
proved by showing that they do not increase the WL-TWSC. The WL-TWSC can be shown to be written as 

WL-TWSC = plht + 2pkhlsl,AkSk,a - 2pkhk^a^ + Pk (19) 

where Ak is the interference plus noise correlation matrix, and Pk is the sum of the correlations not involving 
the signature Sk^a- After replacing Sk^a with dk^a/\\dk,a\\ = — ^ _2 '''''' — -, the WL-TWSC of the new set of 
signatures can be expressed as 

WL-TWSC = 2pkh f'^;ff''^^ +Pk- 2pkhla^ + plkj . (20) 

Next, we compute the EVD of Ak = VAV^ and set x = Sk^a- Then, 

WL-TWSC < WL-TWSC x^A'^x < {x^ A-'^x){x^ Ax) . (21) 



June 18, 2009 



DRAFT 



Exploiting the results reported in the appendix of [4] it can be shown that 

1 < {x^ K-^x){x^ Kx) , and {x^A-^xf < [x^A-^x] . (22) 

Applying these two bounds in succession it can be thus shown that the spreading code update decreases the 
WL-TWSC and, consequently, iterations (fTSl ) admit a fixed point. 

B. Orthogonality of the signatures at the fixed point for K < 2N 

Generalizing the arguments of [4] to the case in which complex spreading codes and WL MMSE filtering 
are adopted, we prove now that for K < 2N, the steady-state output of iterations ([TSll is a set of orthonormal 
augmented signatures, provided that the K augmented signatures used as a starting point of our algorithm are 
linearly independent. Of course, this condition can be met only if the range span of the K starting (non-augmented) 
signatures S{0) = [si, . . . , sk] coincides with C^, with C the complex field. 

Let A be the diagonal matrix whose k-th element on the diagonal is 2pfc/i| = a|, and let S = [si^a, • • • , SK,a]- 
We first prove the following lemma: 

l2N + SaAS^\, Vx>0, (23) 

where Sa is the matrix of the augmented signatures after having updated s,t,a- To prove ( [23l) . let us consider 
first the case x > 0. Set Zj. = J2j^k ^j^j,a.^fa' '^h^'^ 

\xl2N + SaAS^l = \xl2N + Zk + alsk,aS^J = \xI2N + Zk\\l2N + {xl2N + Z k)~-^ k ,aS^J = 

\xI2N + Zk\{l + als^^^{xl2N + Zk)'^Sk,a) ■ 

Similarly we obtain 

\xl2N + SaASa"\ = \xI2N + Zk\{l + a2^<^(xl2iv + (24) 

Substituting for ^ with its expression and taking into account the two previous equations, lemma ( [23] ) can be 
expressed as 

S^^{Zk+2Nol2N)-'iZk+xl2NyHZk+2Nol2N)-'Sk,a ^ H ^ T , '7 ^-l (l')^ 

Upon performing the eigenvalue decomposition (EVD) of = VAV^ and letting w = s^ a, we obtain 

2N I |0 2N : |2 2N , ,0 

^(A, + x)(A, + 2iVo)2 -^A, + x^(A, + 2iVo)2 ^ ^ 
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Since \\w\\ = 1, we are allowed to define a discrete random variable Y with the following probability mass 
function: 



pviy) 

Therefore, ( |26l ) is equivalently written as 



y = h,l<i<N 
otherwise . 



(27) 



> E 



1 



{Y + x) 



E 



1 



{Y + 2iVo)2 



(28) 



For the proof of (1281) refer to [4]. This proves the result for x > 0. In order to prove the result for x = 0, consider 
the following A^th-order polynomial: 



f{x) = \xl2N + SaAS, 



xl2N + SaAS^l = Ul=i{x + A.) - llt=iix + A 



2N 



(29) 



wherein Aj and Aj are the eigenvalues of SaASa and SaAS^ respectively. We have proved that /(x) > 
0,Vx > 0. Then, due to the continuity of f{x), we have /(O) > 0, which proves the lemma for x = 0. 
Given the above lemma, it is straightforward to note that 



H 



xlK + ASa'Sa\>\xlK + AS^Sa\, Vx > . 



Setting X = we obtain 



ASa^'Sal > \AS^Sa 



(30) 



(31) 



Since A is an invertible positive definite matrix, the rank of AS^ Sa equals the rank of Sa- Therefore, if 
the initial matrix Sa has full rank, the matrix Sa at the fixed point will be invertible. 

Let us now prove that at the fixed point the augmented signatures are eigenvectors of the correlation matrix 
of the augmented received vector. At iteration n we can write 

M-\n-l)sk,a{n-l) 



\\M-\n-l)sk,a{n-l)\\ 
At the fixed point everything is independent of n. Therefore 



(32) 



(33) 



and finally we obtain MaS^^a = '^kSk,a^ which proves the result. 

Now, we exploit the above result in order to prove that at the fixed point we obtain an orthonormal set of 
signatures: 

MaSk,a = {SaAS^ + 2Nol2N)Sk,a = <TkSk,a , V/c . (34) 



June 18, 2009 



DRAFT 



10 



Multiplying both sides of (l34l ) by sf^^ we have 

^l^ai^aASa +2Nol2N)Sk,a = CrkS^a^k,a, V/c, L (35) 

Rewriting (1351 ) in a matrix form we obtain 

(5„A5f + 2Nol2N)Sa = S^Sa^ . (36) 

After some algebra we have 

S^Sa = A-\'S-2NolK). (37) 



Since the matrix on the right-hand-side of (1371) is diagonal, at the fixed point we achieve orthogonal augmented 
signatures. Moreover, since the signatures are constrained to have unit norms, we conclude that the algorithm 
converges to an orthonormal augmented set. This is a very attractive feature, since transceiver optimization 
algorithms using linear receivers can provide an orthonormal set of signatures only for a number of users K < N, 
whereas, if WL processing can be used, orthonormal signatures can be granted to a number of users K < 2N: 
this seemingly inexact results is due to the fact that, when WL processing is adopted what really matters is not 
the scalar product among the spreading codes, but its real part. 

1) Optimality of the steady-state solution for k < 2N: Since the signatures are constrained to have unit norm, 
the following bound for the WL-TWSC holds: 

1 ^ 1 

WL-TWSC >^Y."'k = l^^i^^) ■ (38) 

k=l 

It is easily verified that (l38l ) can be achieved when the signatures are orthonormal. Hence, when K < 2N our 
algorithm achieves the global minimum of the WL-TWSC. We also note that optimization of the WL-TWSC 
is equivalent to optimization of other performance measures such as the Total MMSE (TMMSE) [5], and the 
Sum-Capacity (Cgum) [26]. Indeed, the WL-TWSC, TMMSE and Csum for the considered system can be easily 
expressed in terms of the eigenvalues of SaAS^ , [Ai . . . A2Ar] as follows: 



WL-TWSC = 15^ A? , 

1=1 

2N 

C,^ra = Y,\og{l + Xi/2NQ) and (39) 



i=l 

2N , 

TMMSE = K-^— ' 

1=1 



Xi + 2No 

Since SaAS" is an hermitian matrix, its eigenvalues will be real numbers. As a consequence, the WL-TWSC 
and TMMSE are Schur-convex functions, while Csum is a Schur-concave one. Therefore, we conclude that the 
signature set provided by the fixed point of dTSl ) also minimizes the TMMSE and maximizes Csum- 
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C. Properties of the steady-state spreading codes for K > 2N 

So far, we have proved that the algorithm ([TSl l always converges, and that it achieves the global minimum of 
the WL-TWSC for K < 2N. Now, properties of the equilibrium spreading-codes for K > 2N are investigated. 

We assume that the algorithm is initialized with a full rank matrix Sa, so that at the fixed point the matrix 
SaAS^ is positive definite. We start with the following theorem that characterizes the fixed points of ( fTSl ). 

Theorem 1: Let Sa be a fixed point of ( fT5] ). and qi,...,qjy be an orthonormal basis of eigenvectors of 
SaAS^ ■ Then, there exist a number L G {1, . . . , N}, a partition { Ji, . . . , Jl} of the set {1, . . . , K}, a partition 
{/i, . . . of the set {1, ... , N}, and positive numbers Ai, . . . , Al, such that 

X{SaAS^) = [Ai,.^.,Ai^, . . . , Al,.^.,AlJ , (40) 

^^ = ^E"fc' V^G1,...,L, (41) 

Sk,a^Sk,a2 = , Vfci G Je,k2 i Jl and (42) 

{sk,a ■■ k £ Ji} C span{q„ : n £ Ii} , (43) 
where denotes the cardinality of Ii, and the sets J/ and // are defined as follows: 

Ji = {k£l,...,K:SaAS^Sk,a = )^iSk,a} and (44) 

Il = {nGl,...,N : SaAS'^q^ = Xiq^} . (45) 

Proof: A complete proof for the case of real channels can be found in [6]. Here we extend it to the complex 
case, putting emphasis on some points that will be useful in the sequel. Let L be the number of distinct eigenvalues 
of SaASa ■ From previous sections we know that all of the augmented vectors Sfc ^ are eigenvectors of SaAS^ ■ 
By grouping the augmented signatures associated to the same eigenvalues we obtain (l44l ). S ince SaAS„ is an 
hermitian matrix, eigenvectors associated with distinct eigenvalues are orthogonal, thus implying (l42l) . Now, let 
qi,. . . ,q]\f be an orthonormal basis of eigenvectors of the hermitian positive definite matrix SaAS^- If we 
group together the eigenvectors q„ associated with the same eigenvalue, we obtain ( |45l ). The cardinality of 
equals the (geometric and algebraic) multiplicity of the eigenvalue A^. Then, (l40l) is proved. The interested reader 
can refer to [6] for the rest of the proof. ■ 
Summing up, this theorem states that at the fixed point, the matrix SaAS^ has L distinct eigenvalues with 
multiplicities . . . , |/l| respectively, and the columns of the matrix Sa are partitioned into orthogonal sets, 
with each set containing the columns that are eigenvectors associated with the same eigenvalue of SaASa ■ Note 
that since we are assuming a full-rank matrix Sa, and K > 2N, the columns of Sa will span C^^, implying 
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that for each eigenvalue there will be at least a column of Sa that is an eigenvector associated with that 
eigenvalue. This means that we can always choose the orthonormal basis Qi, . . . , so that each contains at 
least one proper augmented vector. Moreover, since K > 2N, and SaAS^ is not singular, denoting by Sj^ the 
matrix whose columns are the vectors belonging to Ji, we have = rank(S'jJ. 

1 ) The suboptimal fixed points of iterations di5D are unstable: Reference [6] proves that for the case of real 
signals and linear processing all the suboptimal fixed points of the MMSE update algorithm are unstable. We 
prove here that a similar result also holds true for iterations ( fTSl ). As ( fTSl ) does not increase the WL-TWSC, all 
we need to prove is that among the fixed points of (fTSl) there are no local minima of the WL-TWSC, except the 
global minimum. 



To begin with, let V and S denote the sets {va : Va 



1,17 G C~}, and {[s 



Sk,a £ V}, respectively. We can define a metric on S as 

d ■■ {S'^, S'a) eS^ ^ d{S'^, S'^) = rnax arccos(s'J„4 ) . 

k=l,...,K 

Note that for vectors in C^^, s'^^s'^^ is in general a complex quantity, and hence d is not a metric. However, 
we now prove that as long as vectors in V are considered, s'^^s'^ ^ is real and d is a metric on S. 

Vi 



Let vi^a 



and V2n 



V2 



be two generic vectors belonging to V. Then we have 



(46) 



which is of course real. Now, if we express vi and V2 in terms of their real and imaginary parts, we obtain 
where vji^i, and i?/ i are the real and imaginary part of the vector i^i; i>i?,2> and i?/ 2 are the real and imaginary part 



of the vector V2 ; vrj 1 = \/2 



vra 



Via 



, and vri^2 = V2 



VR,2 



VT2 



. Therefore, when vectors in V are considered. 



we can always express the hermitian product v^^V2^a in terms of a scalar product between real vectors. Moreover, 
since vi^a and V2^a have unit norm, the vectors vrj^i and vrj^2 will also have unit norm. As a consequence we 
have 

arccos(t;{[^i>2,a) = avccos{v'^j ^^v rj ^2) = avccos{\\v rj ^i\\\\v rj ^2\\ cos(a)) = a G [0,7r], (48) 
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where a is the angle between the vectors f/j/,! and vjii^2- This allows us to prove that d satisfies the triangle 
inequality. In fact, given Sa, S'^, S" G 5 we have 

diSa,S'^) = maxfc=i,...,Karccos(s^„s'' ) = max arccos(s^^ fcS^^ fc) < 

' ' k=l,...,K 



max 
k=l,...,K 



< (49) 



max arccos(Sfc^„Sfc^„) + max arccos(s'fc ,j4',a) = d{Sa, S'^) + d{S'a, S'^) ■ 

ft -L ^ . . . ^ l\ ftl J. J . . . 

The other properties that characterize a metric are straightforward to prove, and hence d is a metric on S. 
In the following we extend to the case at hand some results proved in [6] for real signals. 

Lemma 1: Let G 5 be a fixed point of (fTSl ) and a local minimum of WL-TWSC. Then, given £1,^2 £ 
1, . . . , L with Xi^ > Xg^, ki G J^j, k2 G Je^ we must have a^^ > a^^. 

Proof: Suppose the thesis does not hold. Then, take e > and set a = sin(e), (3 = — ^a. Now consider a 
matrix S'^ with the same columns of Sa except for 



Note that G 5 because s^^ ^ and s^^ q are orthogonal, which implies ^|| = ||s'^^ ^|| = 1. Moreover, a, 
P € TZ. Therefore, Sfc^^a and Sf^^^a are proper augmented vectors. Since we are assuming < a|^, it is easily seen 
that d{S'^, Sa) < e, which means S'^ belongs to the e-neighborhood of Sa- Now we set S'^AS'^ = SaAS^ + A 
where 




(50) 

Then, after some manipulation we get 

Given the hypothesis, and as we are assuming < a|^, for e approaching we have WL-TWSC(5a) > 
WL-TWSC(S^a). Therefore, we have found a matrix in the e-neighborhood of Sa, with a higher WL-TWSC than 
Sa, which contradicts Sa being a local minimum of WL-TWSC. ■ 
Lemma 2: Let S^q G 5 be a fixed point of (fTSl ) and a local minimum of WL-TWSC. Let ^ G 1 , . . . , L with 
Xi > min<./gi_...^i A<?'. Then |J^| < 

Proof: Suppose the thesis does not hold. Then, we can find £1, £2 ^ ■ ■ ■ ,L, with A^^ > A^^ and | J^J > |/^J. 
As a consequence, since |J^J > |/^J = rank(S^j^J, there exists a unit norm vector v G Cl'^^i', belonging to the 
null space of S^j^^ Aj^^ . We now show that, for the case at hand, the null space of Sj^^ Aj^^ also contains a unit 
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norm vector G TZ^'^'^V Indeed, given the structure of the matrix Sj^^Aj^^, and letting M be the null space of 
P = Sj.^Aj^^, we can write 

/ \ / \ 



Then, 



Pi 



E 

1=1 



ViPi = 



(51) 



(52) 



. 



1=1 



Conjugating the second equation, and summing and subtracting, we obtain 



i=l 



(53) 



\, i=i 



Therefore, the vector Vr = 3f?(v)/||3f?(i>)||, is a real unit norm vector belonging to A^, whose components we 
denote by [vr,i, ■ ■ ■ ,Vr^k, ■ ■ ■ !^r,|jfj|]- Now take e > and define a matrix S'^ with the same columns of Sa 
for k ^ and s'j, ^ = (cosafc)sfc^a + (sinafc)q, for k G J^^, where = ew^^fc, and q £ VPl-^^s- pointed 
out in the previous section, such a vector always exists. Therefore, s'^ a ^ which implies S'^ G S. Note also 
that d{Sa, S'^) = emaxfcgj^^ \vr,k\ < e> implying that S'^ is in the e-neighborhood of Sa- Again, perturbing the 
signature matrix, we obtain WL-TWSC(5a) - WL-TWSC(S';) = i {2e^{Xi, - Xi,)\\Aj^^Vrf + o{e^)). Since 
for e approaching we have WL-TWSC(Sa) > WL-TWSC(5^), we can find a matrix in the e-neighborhood of 
Sa with a larger WL-TWSC, which contradict Sa being a local minimum of WL-TWSC. ■ 
Similarly to [6], these results can be exploited to prove that among the fixed points of (fTSl) there are no 
local minima of the WL-TWSC except the global minimum, implying that all the non-optimal fixed points of our 
algorithm are unstable. Hence, a noisy version of (fTSl ). obtained by adding a small perturbation to the deterministic 
updates, will converge to the global minimum of the WL-TWSC with probability 1. 



D. Oversized users and WL filtering 

Reference [2] shows that, for the case in which K > N and real signatures with linear detection are considered, 
the optimal spreading code allocation may grant to some large-power users orthogonal channels. Otherwise stated, 
if one or more users are received with a power that is large relative to that of the remaining ones, these users. 
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called oversized users, are given orthogonal spreading codes and enjoy single-user performance; the remaining 
users, instead are given generalized Welch-Bound-Equality (WBE) sequences and are confined in the orthogonal 
complement of the subspace spanned by these oversized users. More precisely, in [2], user i is defined to be 
oversized if 



dl > .'r'~J.K ^y^' (54) 

where df = pihf \/i = I, . . . , K, and ldi>di^ = 1 if > and zero otherwise. Obviously, there can be at 
most — 1 oversized users. 

Now, let us focus on the WL filtering case with complex spreading codes. We have seen that for K < 2N 
it is possible to grant each user a spreading code that is mapped into an orthogonal augmented signature at the 
receiver, which is of course the optimum spreading code allocation. As a consequence, for < K < 2N the 
distinction between oversized and non-oversized users is meaningless if WL filtering is used, while it is necessary 
in the linear filtering case. Consider now the more interesting case in which K > 2N. When WL filtering is 
considered, it can be shown that a user is oversized if 

> ' i ■ (55) 



- Ef=l la?>af 



Note that now there can be at most 2N — 1 oversized users. Again, the optimum spreading code allocation is 
to grant each oversized user an orthogonal augmented spreading code, while other users are given generalized 
WBE sequences. 

Now, comparing relations (l54l ) and dSSl ) it can be shown that the set of oversized users resulting from iterations 
([TS] ) includes the set of oversized users resulting from the use of iterations ( fT6l ). Indeed, assuming user i satisfies 
(l54l ). we can write 

,2 YL%\d]'^d\>d] ,2 ^%\d]\d\>d] 
df > — - — — => df > ^ — ^ 

N - Ef=i ld?>d? 2N - Ef=i ^>d^ 

K K 2 (^^^ 

^ 2df > — - — — ^ af > ^ — — . 

Accordingly, denoting by /Cl and /CwL the sets of oversized users in the linear and WL case respectively, we 
have the relation /Cl ^ /CwL- Of course, having more oversized users is a good occurrence, since more users 
will have orthogonal channels. 

E. Equivalence between the minimum-MSE and the maximum-SINR non-cooperative games 

Equation (ITll) shows the link between the fc-th user's MSE and its achieved SINR under the assumption that 
the WL MMSE receiver is used. Given this one-to-one strictly decreasing relationship, it thus follows that non- 
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cooperative MSE minimization is equivalent to non-cooperative SINR maximization. A different, yet equivalent 
form for the game (fT3]) is thus 

max 7fc , k = 1, . . . , K — 1 . (57) 



F. Sum-capacity of a CDMA system using as spreading codes the fixed point of iterations di5D 

The sum-capacity of a CDMA system with processing gain N using the spreading codes resulting from the 
unique stable fixed point of iterations (flSl ) can be contrasted to that of a CDMA system using complex 2N- 
dimensional codes with linear detection. In this section we thus compare the sum-capacity of the following two 
systems: 

a) a CDMA system with processing gain and a WL receiver. This system can be seen as a CDMA system 
employing 2A^-dimensional, unit-norm, augmented spreading codes, that is codes in S. 

b) a CDMA system employing 2A^-dimensional, unit-norm, complex spreading codes and linear detection. 
For the sake of comparison, we suppose that the noise power, the channel gains, and transmit powers are the same. 
Therefore, the parameter a, and the matrix A are the same for both systems. Now, denote by Sc and c\^^^ 
the spreading matrix of system b), and the set of all matrixes in c2AfxX unit-norm columns, respectively. 
The sum-capacities of the two systems can be written as [1] 

Csum,WL = max log|/2Af + cr"^5aA5f I and 

(58) 

Csum,c = max log|J2Af + o- ^ScAS'f I . 

Since S C Cf^^^ , it is clear that Csum,WL < C'sum,c- However, we will now prove that Csum,WL = C'sum,c> 
and therefore no loss in performance (in terms of sum-capacity) is suffered when using codes in S. To see this, 
we will consider separately the cases K < 2N and K > 2N. Let us first assume that K < 2N. We can write 

Csum.c = ^ max log|l2Ar + a''^ScAS^\ = max log|Ii^ + a~'^S^ScA\ = log\lK + cJ^'^A] (59) 

In the last equality we have exploited the result that for K < 2N the maximum is achieved by orthonormal 
signatures. As for Csum.WL we have 

Csum,WL = max log|/2Ar + (J''^SaAS"\ = max logllic + a~'^S^SaA\ . (60) 
Saes Saes 

Again, the maximum is achieved by a matrix with orthonormal columns. From previous sections we know that 
such a matrix is always to be found in S for K < 2N, and therefore we can write 

Csum,WL = max log|/i^ + fJ^^AI = Csum,C • (61) 

Saes 
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Consider now the case K > 2N . First, we express the sum-capacity of both systems in terms of the vectors 
X{ScAS^) and \{SaAS^), whose components are the eigenvalues of the matrix ScAS^ and SaAS^, 
respectively: 

2N 

" " (62) 



Csum,c = max y^log(l + (Ai/cr^)) and 



2N 



Csum.WL = max y'log(l + (Aj/o-^)) . 



(63) 



In order to proceed with the proof, we need a preliminary result first. We need to prove that Csum.WL = 2Csum,R^ 
where Csum,R is the sum-capacity of a CDMA system employing 2A'^-dimensional, real, unit-norm, spreading 



codes. Denoting by Sr and TZ\^^^ the real-valued spreading matrix, and the set of all matrixes in IZ 



2NxK 



with 



unit-norm columns, respectively, we have 

1 



Csum,R — „ max — 



log\l2N + (r-'^SrASj\ 



max — 



log|/j^ + a-^Sj.SrA\ . 



(64) 



Note that [5], for real spreading matrixes the sum-capacity is defined with a factor 1/2. 



Now, consider the following bi-injective transformation between the sets S and TZ^ 



2NxK 



f:Sa 



VI...VK 



vi . . .V 



K 



^{vi)...^{vk) 



2NxK 



(65) 



This implies that for any matrix Sa G S, there exists a matrix Sr = f{Sa) & TZ\^^^ such that 
Therefore we can write 

max log\lK + cr~'^SrSrA\ = 2Csum,R ■ 



Csum.WL = max log|Iii' + fJ "^S^SaAl = 

5„g5 S^enl"^" 



(66) 



(67) 



Bearing this result in mind we can resume our main proof. In fact, generalizing the arguments of [2] to the case 
of complex 2A^-dimensional spreading codes, it can be shown that both the eigenvalues \{ScAS^) as Sc takes 
value in Cl^^^ , and the eigenvalues X{SrAS^) as Sr takes value in 1Z\^^^ , define the set 

{(Ai, . . . , A2Ar) G 7^^^ : (Ai, . . . , A2Ar,0, . . . ,0) majorizes (oi, . . . , Oi^)} , 

which impUes Csum.c = 2Csum,R- Thus, by virtue of ([67]), we get Csu^.c = 'iCsnm,K = Csum,WL- 

Therefore, a CDMA system with processing gain N, a WL receiver and spreading codes resulting from 
iterations (fTSl ). has the same sum-capacity of a CDMA system with processing gain 2N, a linear receiver and 
spreading codes resulting from iterations ([T6l ). 
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IV. Non-cooperative power and transceiver optimization for maximum energy-efficiency 

So far, we have considered maximum-SEMR and minimum-MSE games with respect to the spreading code and 
WL receiver of every user. Now, we include the transmit power among the parameters to be tuned, and consider 
the non-cooperative power control and WL transceiver design game assuming that the objective function is 
the energy efficiency, i.e. the number of bits reliably delivered to the receiver for each energy unit used for 
transmission. 

To be more definite, assume that each mobile terminal sends its data in packets of M bits, and that it is 
interested both in having its data received with as small as possible error probability at the AP, and in making 
careful use of the energy stored in its battery. Obviously, these are conflicting goals, since error-free reception 
may be achieved by increasing the transmit power, which of course comes at the expense of battery life. A useful 
approach to quantify these conflicting goals is to define the utility of the fc-th user as the ratio of its throughput, 
defined as the number of information bits that are received with no error in unit time, to its transmit power [16], 
[17], i.e. 

Uk = Tk/pk ■ (68) 

Note that is measured in bits/Joule, i.e. it represents the number of successful bit transmissions that can be 
made for each battery energy-unit used for transmission. The utility function (1681 ) is widely accepted and indeed 
it has been already used in a number of previous studies such as [15]-[18]. Denoting by R the common rate of 
the network and assuming that each packet of M symbols contains L information symbols and M — L overhead 
symbols, reserved, e.g., for channel estimation and/or parity checks, and following the reasoning of [15], [16], a 
faithful and mathematically tractable approximation for the utility Uk in ( [68l ) is the following: 

Uk = R^^^, yk = l,...,K. (69) 

In the above equation, /(7fc) is the so-called efficiency function, approximating the probability of successful 
(i.e. error-free) packet reception. As an example, for binary phase-shift keying (BPSK) modulation, the choice 
filk) = (1 ~ e""^*")*^ is widely accepted. The results of this paper, however, hold not only for this particular 
choice, but for any efficiency function /(•) that is increasing, S-shaped, approaching unity as —>■ +oo, and 
such that /(7fc) = 0(7^) for vanishing 7^. 

We consider now several non-cooperative games for utility maximization. First of all, we note that the games 
reported in [15] (optimization with respect to transmit power only), in [16] (optimization with respect to transmit 
power and linear receiver choice), and in [17] (optimization with respect to transmit power, hnear receiver choice 
and spreading code choice) can be extended to the case that a complex signal, as the one in ([T]), is considered. 
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We do not give here the full details to avoid duplication of these previous results, but simply mention that it can 
be shown that for all of these games a unique NE point exists, and that, for the case in which also spreading 
code optimization is performed, the NE point is Pareto-optimal for K < N. We prefer to give much more details 
on the interesting case in which WL filtering is used. 

Let us thus assume that the decision rule is adopted, and consider the problem of utility maximization with 
respect to the transmit power and the WL receiver choice, i.e. 



with pk < Pk. max the maximum allowed transmit power for the k-th user. We have the following result. 
Proposition 1: The non-cooperative game defined in d70D admits a unique NE point (p^, d^), for k = 1, . . . ,K, 
wherein 

- d^. is the unique vector, up to a positive multiplicative constant, resulting from the upper N -entries of the 
k-th user widely linear receive filter; and 

- Pk = mm{pk, -Pfc.max}. with pk the k-th user's transmit power such that the k-th user's maximum SINR 7^ 
equals 7, i.e. the unique solution of the equation 7(7) = 7/' (7). with /'(7) the derivative of f{'~f). 

Proof: The proof can be given following the same arguments as in [16]. The full details are however omitted 
for the sake of brevity. ■ 
Consider finally the case in which also spreading code optimization is performed, i.e. we have 



We now have the following result. 

Proposition 2: The non-cooperative game defined in ([77]) admits a unique NE point {p^, d^, s^), for k = 1, . . . , K, 

wherein 

- s|. and d|. are the unique (up to a positive scaling factor for the linear receiver; unicity for the spreading 
codes set is instead meant with respect to their correlation matrix) k-th user's spreading code and receive 
filter resulting from the upper N -entries of the fixed points of iterations 



max Ufc , VA; = 1, . . . , , 

Pk,dk 



(70) 



max Uk 
Pk,dk,Sk 



\/k = l,...,K . 



(71) 




,a 




K 



< 




k = l,...,K 



Denote by 7^ the corresponding SINR. 
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- p'f. = min{pfc, Pfc,max}> with pk the k-th user's transmit power such that the k-th user's maximum SINK 7^ 
equals 7, i.e. the unique solution of the equation /{"f) = 7/' (7), with f'{"f) the derivative of /{"f). 

Moreover, for K < 2N, the NE point is Pareto-optimal. 

Proof: Let D = [di^a, ■ ■ ■ ,dK,a], and X = {S,D), and denote by = [pi, . . . ,Pk-i,Pk+i, ■ ■ ■ ,PkV the 
{K — 1) -dimensional vector containing the transmit powers of all the users except the k-th one. Upon defining 

1 



i^k 



(72) 



" 2hl{dl,Sk,a? 

the utility function for the considered game can be written in the form 

Uk{pk,Ik{X,p_k)) . (73) 

Since, for fixed power p^, the utility function is a decreasing function of I^, the utility function ( |69l ) is said to be 
separable in the two parameters, p and X [27], and the corresponding game is a separable gamea Let us denote 
by Qx{p) and by Qp{X) the subgame arising from utility maximization with respect to the spreading code and 
uplink linear receiver optimization for a fixed transmit power configuration p, and the subgame arising from utility 
maximization with respect to transmit power for fixed spreading codes and linear receivers, respectively. Based on 
the results of the previous section it can be shown that the subgame Qx{p*) admits a NE point X* = {S* , D*), 
arising from the fixed point of the iterations ([TSl l. Similarly, for any user k, it can be shown that the transmit 
power p^,, as defined in the text of Proposition 2, is an NE point for the subgame Qp{S* , D*) [15]. Given any 
interference Ij^{X ,p_i^), the power = min(7/fc(X,p_fc), Pjt max) maximizes the utility of the user k and it 
is a continuous function of the interference Ik. Moreover it is easily seen that denoting by X*{p) the NE of 
the subgame Qx{p), the interference function Ii^{X*{p),p_i^) in ( 1721 ) is continuous in p, for k = 1, . . . , K. 
According thus to Theorem 2 in [27] the existence of an NE for the game (TtT] ) is guaranteed. Finally, as regards 
the fact that the NE is Pareto-optimal for K < 2N, this immediately descends from the result, shown in the 
previous section, that for K < 2N iterations ( [T5l ) converge to a set of orthonormal augmented signatures, thus 
implying that the users at the NE enjoy separate channels, with no multiuser interference. ■ 

A. Network performance prediction through LSA 

Following the approach of [17], we now show how LSA arguments can be used to predict the network 
performance. In particular, we focus here on predicting the utility profile of the active users in a large CDMA 

^Otherwise stated, in a separable game the utility function can be written as in il3i . and, for any fixed pk, the utility function is a 
decreasing function of It. 
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system (i.e. a CDMA system with N and K — > oo but K/N = a) employing WL filtering at the receiver, 
and with no spreading code optimization. LSA tools were first developed in [28], where the following heuristic 
relation for the SINR achieved by the generic k-th user with linear MMSE detection and random unit-norm 
spreading codes is derived for real channels 

Generalizing (1741 ) to the case of a complex channel and WL filtering at the receiver, the SINR 7^ achieved by 
user k can be shown (the proof is omitted for the sake of brevity) to satisfy the equation 

^ '^Pkhl 

r) AT I 1 ^Vkll^PjIlj 

The above equation can be used to come up with a non-iterative power control equation for the maximization of 
the energy-efficiency. Indeed, since all users are to achieve the same target SINR 7, it is reasonable to assume 
that they must be received with the same power, i.e. 

2hlpi = 2hlp2 = ... = 2hlpK = Pr . (76) 

Substituting (1761 ) into (1751 ) we obtain 

Ph^^^. (77, 



2(1+7) 

Therefore we can devise the following power control algorithm 

Pfc = min [ ,Pmax ) V/c = 1, . . . , (78) 

with Pmax the maximum allowed transmit power, assumed to be the same for all the users. Of course, once the 
transmit powers are known, the actual achieved SINRs can be computed using (ITSl ). and the achieved utilities 
are given by ( l69l ). Equation (1781 ) is the generalization (to the case of WL filtering) of a power control algorithm 
reported in [28]. However, it assumes that all users are received with the same power, which is a realistic 
assumption only if no user ends up transmitting at maximum power. In fact, a user transmitting at maximum 
power will not achieve the target SINR and therefore (1761 ) is no longer true. Consequently, as our numerical 
results will show, ( TTSl ) exhibits a poor prediction accuracy when many users are transmitting at their maximum 
power. To circumvent this problem, note that in order to carefully predict the users' utility profile, the number 
of users transmitting at maximum power needs to be estimated. Denoting by Nm such a number, an estimate of 
Nm can be obtained according to the following formula: 

^ = E ^ [^-r^%r- - Pmax) , (79) 
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with ti(-) denoting the unit-step function. It is also fairly reasonable to assume that the users transmitting at 
maximum power will be the Nm users with the smallest modula of the channel coefficients. Therefore, ordering 
the modula of the channel gains in descending order we can write now the SINR for the generic active user as 
follows: 

2NPk 

ANN -U ^^P" -U 2P,p^,^hl - T ' ^^^^ 

with Pk = 2hf.pk. Solving this equation with respect to gives the desired receive power for user k. Actually, 
since equation (l80l) can be written for any k = I, . . . ,K, the solution is the desired receive power for all 
users, that is Pr. Once Pr is known, the powers for all users are easily found as follows: 

Pk = min(Pfc/2/ifc,Pmax) , k = I, . . . ,K . (81) 

Once the powers are known, the actual achieved SINRs, can be easily computed for all users. Then, the users' 
utilities are given by 

'"k = irf > k=l,...,K. (82) 

V. Numerical results 

In this section we present some simulation results that give insight into the performance of the proposed 
transceiver optimization algorithm and, also, corroborate, the validity of the theoretical findings. We consider 
a DS/CDMA system with BPSK modulation (so that the received signal is improper), and with a randomly 
generated starting signature set. 

In Figs. 1 and 2 we compare the transceiver optimization iterations ([TSl l with the classical MMSE signature 
update iteration of [4]. The processing gain has been set to iV = 15, and the users' received powers have been 
randomly generated. Denote by Sa, and S the sets of augmented and non-augmented signatures, respectively. 
Fig. 1 shows the minimum and maximum eigenvalue of the matrices Sa, and S versus the iteration index 
for two different number of users, namely K = 10 and K = 20. As expected, for K = 10, the maximum and 
minimum eigenvalues converge to 1 in both cases, implying that both S, and Sa converge to the identity 
matrix Ik, whereas for K = 20 only Sa Sa converges to Ik- This can be explained noting that since the matrix 
S^ S has rank A^, it cannot converge to Ik when K > N, whereas the matrix Sa Sa has rank 2A^, which 
makes the convergence to Ik possible also for A^ < < 2A^. Fig. 2 shows the WL-TWSC and the TWSC 
of the non-augmented set S versus the iteration index for K = W and K = 20. Both the WL-TWSC and the 
TWSC can be lower bounded as in ( [38] ). Given the randomly generated received powers the lower bounds have 
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been found to be 

tr(A2) = 5.36 for K = 10 

(83) 

tr(A2) = 12.08 for K = 20 

As expected, the lower bound is achieved by both WL-TWSC and TWSC when K = 10, but only by WL- 
TWSC when K = 20. These two figures give evidence of the fact that WL filtering coupled with spreading 
code optimization can double the number of users able to enter the network without suffering any multiuser 
interference, with respect to previous transceiver optimization algorithms. 

Fig. 3 addresses the performance of iterations ([TSl l when K > 2N. Here, the processing gain has been set to 
N = 5, the number of users to K = 12, and the users' received powers are 

Pi = 11.51 , P2 = 7.94 , and P^ = 1 for k = 3, . . . , K . (84) 

It is easy to see that the two users with the largest powers are oversized [2]. In this scenario, the optimum 
signature set dedicates each oversized user a signature orthogonal to all other users, while the other users are 
given generalized Welch-bound signatures. Given the powers in (l84l ). the eigenvalues of SaAS^ associated to 
the optimum set of signatures can be shown to be 

P 

Ai = Pi, X2 = P2, and Afc = |^^^ = 1.25 for fe = 3,...,K. (85) 

Fig. 3 shows the eigenvalues of SaAS^ versus the iteration index. As expected, the eigenvalues converge to the 
optimum values, implying that Sa converges to the optimum set. 

Consider now the system performance of the maximum energy-efficiency non-cooperative game at the NE. 
Assume a processing gain = 11, and a packet length M = 120; for this value of M the equation /{j) = 7/'(7) 
can be shown to admit the solution 7 = 6.689 = 8.25dB. A single-cell system is considered, wherein users may 
have random positions with a distance from the AP ranging from 10m to 500m. The channel coefficient h^e^'f"' 
for the generic k-th user is assumed to be a zero-mean complex Gaussian random variate with variance equal to 
, with dk being the distance of user k from the access point (AP). We take the ambient noise level to be 
Mo = 2.5- 10"^°W/Hz, while the maximum allowed power Pfc^max is OdBW. We present the results of averaging 
over 10^ independent realizations for the users locations, fading channel coefficients and starting set of spreading 
codes. More precisely, for each iteration we randomly generate an x A'-dimensional spreading code matrix 
with entries in the set | — 1/ 1 / | ; this matrix is then used as the starting point for the games that include 
spreading code optimization, and as the spreading code matrix for the games that do not perform spreading code 
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optimization. Note that even though the starting spreading codes are real, at the fixed point of iterations (fTSl) 
they are complex vectors. 

Figs. 4-7 show the achieved average utility (measured in bits/Joule), the average user transmit power, the 
average achieved SINR, and the average fraction of users transmitting at maximum power at the receiver output 
versus the number of users, for the game in [16], the game in [15], the game in [17] and for the newly proposed 
non-cooperative games with WL processing at the receiver. Inspecting the curves, it is seen that the approach 
based on WL filtering largely outperforms the games of [15]-[17]. In particular, it is seen that for K < 2N a very 
substantial performance gain can be obtained by resorting to spreading code optimization and WL processing; 
indeed, when K < 2N, the use of spreading code optimization coupled with WL filtering permits granting to 
the active users separate (i.e., orthogonal) channels, so that the multiaccess channel reduces to a superposition 
of K separate single-user additive white Gaussian noise (AWGN) channels. Evidence of this is given by the 
horizontal lines representing the system performance of spreading code optimization plus WL filtering in the 
range 1 < K < 2N. In Fig. 7 we show the fraction of users transmitting at the maximum power: this is actually 
the fraction of users not being able at the NE to achieve the target SINR. As expected, the smaller fraction 
corresponds to the games based on the WL filtering, and also in this case we have an horizontal line representing 
the performance for K < 2N. 

Figs. 8 compares the performance of the LSA prediction algorithms (178] ) and dSTT ). The actual users' utility 
profiles obtained through power control and both WL and linear receiver optimization have been contrasted to the 
users' utility profiles predicted by (178] ) and (ISTT ). It is seen that both with WL and linear filtering at the receiver, 
the LSA-based algorithm is capable to predict the actual users' utility profile with a much higher accuracy than 
the algorithm (iTSl ). In particular, note how the prediction accuracy of (iTSl ) gets worse as the number of users 
increases, while the accuracy of the newly proposed algorithm is almost insensitive to the number of users. This 
can be explained noticing that a larger fraction of users ends up transmitting at maximum power as K increases, 
and recalling from section IIV-AI that (iTSl ) does not take this circumstance into account. 

VI. Conclusions 

This paper has considered the issue of non-cooperative transceiver optimization in wireless data networks 
with WL filtering. Unlike previous work in this area, a baseband complex representation of the data has been 
considered, and emphasis has been given to the case in which the received signal is improper, showing that 
the system performance can be significantly improved by resorting to WL filtering structures. We have solved 
the problem of non-cooperative spreading code design and WL receiver choice for SINR maximization and/or 
MSE minimization, and have provided an update algorithm (i.e., iterations ( [TS] )) that has one stable fixed point 
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coincident with the global optimum point. Relevant properties of the fixed point of Eq. (fTSl ) have been shown, 
discussed, and confirmed by numerical simulations. In particular, spreading code optimization coupled with WL 
filtering has been shown to allow complete suppression of multiuser interference even in overloaded networks, 
up to a number of users that is twice the processing gain. It has been also shown that WL filtering brings 
a substantial improvement to the energy-efficiency of multiuser wireless data networks, and the corresponding 
system performance at the NE has been analyzed, also through LSA arguments. 
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Fig. 1. Maximum and minimum eigenvalues of Sa and S S versus the iteration index for K = 10, 20. The system processing gain 
is TV = 15. 
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Fig. 2. WL-TWSC and TWSC versus the iteration index for K = 10, 20. The system processing gain is A'^ = 15. 
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N=5,K-12. Eigenvalues of SgAS^; 
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Fig. 3. Eigenvalues of SaASa versus the iteration index for K = 12. The system processing gain is TV = 5. 
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Fig. 4. Achieved average utihty versus number of active users for the maximum energy-efficiency non-cooperative games. The system 
processing gain is TV = 11. 
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Fig. 5. Average transmit power versus number of active users for the maximum energy-efficiency non-cooperative games. The system 
processing gain is TV = 11. 
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Fig. 6. Achieved average output SINR versus number of active users for the maximum energy-efficiency non-cooperative games. The 
system processing gain is = 11. 
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Fraction of Users Transmitting at IVlaximum Power. N-1 1 . Power decay factor n-3 
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Fig. 7. Average fraction of users transmitting at tiieir maximum allowed power versus number of active users for the considered 
non-cooperative games. The system processing gain is TV = 11. 
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Fig. 8. LSA-based users' utility profile. The system processing gain is A*' = 64. The power decay factor is n=3 
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